Parametric control of a superconducting flux qubit 



S. Saito,i'2 T. Meno,3 M. Ueda/-^-^ h. Tanaka/-^ K. Semba,!'^ and H. Takayanagi^-^ 

^ NTT Basic Research Laboratories, NTT Corporation, Atsugi, Kanagawa 243-0198, Japan 
^ CREST, Japan Science and Technology Agency, Kawaguchi, Saitama 332-0012, Japan 

^ NTT Advanced Technology Corporation, Atsugi, Kanagawa 24-3-0198, Japan 
'^Department of Physics, Tokyo Institute of Technology, Meguro, Tokyo 152-8551, Japan 

Parametric control of a superconducting flux qubit has been achieved by using two-frequency mi- 
crowave pulses. We have observed Rabi oscillations stemming from parametric transitions between 
the qubit states when the sum of the two microwave frequencies or the difference between them 
matches the qubit Larmor frequency. We have also observed multi-photon Rabi oscillations cor- 
responding to one- to four-photon resonances by applying single-frequency microwave pulses. The 
parametric control demonstrated in this work widens the frequency range of microwaves for control- 
ling the qubit and offers a high quality testing ground for exploring nonlinear quantum phenomena. 
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Quantum state engineering has become one of the 
most important arenas in quantum physics. In particu- 
lar, the coherent control of quantum two-state systems 
(TSS), which are applicable to quantum bits (qubit), 
has attracted increasing interest in the context of quan- 
tum computing and quantum information processing 
Various candidate physical systems are being stud- 
ied for the future implementation of qubits. These in- 
clude artificial quantum TSS like superconducting qubits 
i^lEQtlEIEQilEI^ as well as naturally existing quantum 
TSS like nuclear spins [lol lll| . Of the many candidates 
that may enable us to realize quantum computation, su- 
perconducting qubits based on Josephson junctions have 
gained increasing importance because of their potential 
controllability and scalability. 

The coherent control of a single qubit has been demon- 
strated in many types of superconductiiig circuits, such 
as charge 01, charge-phase 0, phase 0, 0, and flux 
qubits 0, B ■ Recentl y, tw o-qubit operation has been 
demonstrated in charge and phase qubits 0|. In 
addition to qubit operation, the superconducting qubit 
offers a testing ground for exploring interactions be- 
tween photons and artificial macroscopic objects, which 
we shall refer to as "atoms". In the weak-driving limit, 
the interaction between a single "atom" and a single mi- 
crowave photon has been demonstrated with a charge 
qubit, which is strongly coupled to a superconducting 
transmission line resonator |l4| . In the strong-driving 
regime, superconducting qubits have exhibited nonlinear 
optical responses: multi-photon Rabi oscillations have 
been observed in a charge qubit by using microwave 
pulses jl^j l , and under continuous microwave irradiation, 
multi-photon absorption has been observed in a phase 
[l^ and a flux qubit 

In this Letter, we describe the parametric control of 
a superconducting flux qubit with two-frequency mi- 
crowave pulses. We have succeeded in observing two- 
photon Rabi oscillations of the qubit caused by a para- 
metric transition when the qubit Larmor frequency 



matches either the sum of the two microwave frequen- 
cies or the difference between them. We also show 
multi-photon Rabi oscillations corresponding to one- 
to four-photon resonances under single-frequency mi- 
crowave pulse irradiation where the qubit Larmor fre- 
quency was equal to multiples of the microwave fre- 
quency. The parametric and multi-photon transitions 
clearly exhibit high nonlinearity due to interactions be- 
tween a single "atom" and microwaves. In addition, they 
can be useful when designing practical gate operation on 
superconducting qubits. For example, two-photon pro- 
cesses can be used to entangle two qubi ts as already 
shown in the field of ion trap qubits [ijj- Thus, non- 
linear effects will offer new possibilities for flexible qubit 
control as well as being of fundamental physical interest 
in the field of superconducting qubits. 

If we are to observe highly nonlinear phenomena in 
superconducting qubits with strong microwave driving, 
the electromagnetic environment of the qubits must be 
well controlled because the strong driving easily excites 
unwanted environmental resonances, which destroy the 
qubit coherence. Our device was fabricated using elec- 
tron beam lithography and shadow evaporation tech- 
niques defining an inner aluminum loop forming the qubit 
and an outer loop enclosing the dc-SQUID loop used for 
the readout (Fig. ^a)). The inner loop contains three 
Josephson junctions, one with an area a ~ 0.8 times 
smaller than the nominally identical area of the other two 
with critical current Ic ~ 430 nA. The outer loop com- 
prises two Josephson junctions of critical current « 140 
nA. We placed an on-chip microwave line close to the 
qubit at a distance of 20 /im so that the qubit could be 
strongly driven by oscillating magnetic fields, which are 
induced by microwave currents fiowing through the line. 
To control the electromagnetic environment surrounding 
the qubit, we put resistors Rn, i?vi, lead inductances L, 
and shunt capacitors C, 2C on the chip (Fig.H^b)). The 
resistors damp unwanted resonances, which are generated 
in, for example, leads and parasitic capacitors outside the 
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FIG. 1: (a) Scanning electron micrograph of a flux qubit 
(inner loop) and a dc-SQUID (outer loop). The loop sizes 
of the qubit and SQUID are 10.2 x 10.4 p,m^ and 12.6 x 
13.5 /im^, respectively. They are magnetically coupled by 
the mutual inductance M ft; 13 pH. (b) A circuit diagram of 
the flux qubit measurement system. On-chip components are 
shown in the dashed box. L k, 140 pH, C ft; 9.7 pF, _Rii = 
0.9 kfi, _Rvi = 5 kSl. Surface mount resistors i?i2 = 1 kf2 and 
i?v2= 3 kf2 are set in the sample holder. We put adequate 
copper powder filters CP and LC fllters F and attenuators A 
for each line. 



resistors. There are two well-controlled resonance modes 
that are produced in the circuits inside the resistors, 
namely, the on-chip components close to the qubit. One 
is the dc-SQUID's plasma mode with a frequency of 1.0 
GHz, which is formed in the two symmetrical loops, each 
being composed of L, C, 2C, and the SQUID's Joseph- 
son inductance. The other is the harmonic LC resonance 
mode with a frequency of 4.311 GHz, which is produced 
in the larger loop consisting of the two L's and the two 
Cs. In this way, we have achieved an artificial TSS in a 
well-controlled environment. 

The three Josephson junctions of the qubit form a 
double-well potential in the space of the Josephson phase 
when about half a flux quantum threads the qubit loop. 
We use the two lowest levels in the potential as the 
qubit states, which are well separated from the higher 
levels. Thus, the qubit is described by the Hamilto- 
nian iJqb = {Ti/2){eaz -t- AtTa.), where a^.z are the Pauli 
spin matrices. The eigenstates of az describe clock- 
wise and counter-clockwise persistent currents in the 
qubit. The qubit tunnel splitting is described by ?iA, and 
fie = 2/p$o(*i'qb/^'o — 1-5) is the energy imbalance be- 
tween the two potential wells caused by the externally ap- 
plied magnetic flux threading the qubit loop $qb, where 
$0 = h/2e is the flux quantum and /p — /ca/I — {l/2aY 
is the magnitude of the qubit persistent current when the 
qubit is in the tr^ eigenstates and Jc is the critical current 
of the larger junctions. The energy difference between the 
ground state |g) and the first excited state |e) of the qubit 
is ?iWqb = ^/qb = fi^/ A^. Assuming that the applied 
microwaves are in coherent states, we may describe the 
qubit under microwave irradiation by the Hamiltonian 

^ ^ 2 ^'^'^^ ^ + ^ "^^9 kakCTz cos ujMWkt, (1) 
fc=i 



where I is one (two) in the case of a single- (two-) fre- 
quency microwave and is the coupling between the 
qubit and the fc-th microwave (MWfc), whose ampli- 
tude and frequency are and fuwk = ^MWfe/27r, re- 
spectively. Solving the Schrodinger equation with the 
Hamiltonian ^ without the rotating wave approxima- 
tion, we obtain the time evolution of the probability 
Pc(t) with which we find the qubit in |e). The probabil- 
ity Pc (t) oscillates periodically under resonant conditions 
^qb = '^I'j-'MWi (for the single-frequency microwave) or 
^qh = l^^i^Mwi i n2a;MW2| (for the two- frequency mi- 
crowave), where is the MWfc photon number. When 
we operate the qubit away from its degeneracy point 
e ^ 0, the frequency of the oscillation is given for the 
single- frequency microwave irradiation by 
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and for the two-frequency microwave irradiation by 
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Here J„^. is the rifc-th order Bessel function of the first 
kind and A = s/ujqb ~ 1. This approximation is valid 
when e >> A. The dressed atom approach gives results 
similar to Eq. © HI [3 . 

The measurements were carried out in a dilution refrig- 
erator. The sample was mounted in a gold plated cop- 
per box that was thermalized to the base temperature 
of 20 mK {k^T << huiqb)- To produce two-frequency 
microwave pulses, we added two microwaves MWl and 
MW2 with frequencies of /mwi and /mw2, respectively 
by using a splitter SP (Fig.^b)). Then we shaped them 
into microwave pulses through two mixers. We mea- 
sured the amplitude of MWfc 14iwfc at the point between 
the attenuator and the mixer with an oscilloscope. We 
confirmed that unwanted higher-order frequency compo- 
nents in the pulses, for example |/mwi ± /mw2|, 2/mwij 
and 2/mw2 are negligibly small under our experimen- 
tal conditions. First, we choose the operating point by 
setting <I>qb around I.S^&q, which fixes the qubit Larmor 
frequency /qb. The qubit is thermally initialized to be 
in |g) by waiting for 300 /j,s, which is much longer than 
the qubit energy relaxation time (for example 3.8 /is at 
/qb = 11-1 GHz). Then a qubit operation is performed 
by applying a microwave pulse to the qubit. The pulse, 
with an appropriate length tp, amplitudes Vuwk, and 
frequencies fuwk, prepares a qubit in the superposition 
state of |g) and |e). After the operation, we immediately 
apply a dc readout pulse to the dc-SQUID. This dc pulse 
consists of a short (70 ns) initial pulse followed by a long 
(1.5 lis) trailing plateau that has 0.6 times the amplitude 
of the initial part. For $qb < 1.5$0j if the qubit is de- 
tected as being in |e), the SQUID switches to a voltage 
state and an output voltage pulse should be observed; 
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FIG. 2: Experimental results with single-frequency microwave 
pulses, (a) Spectroscopic data of the qubit. Each set of the 
dots represents the resonant frequencies /res caused by the 
one to four-photon absorption processes. The sohd curves are 
numerical fits. The dashed line shows a microwave frequency 
/mwi of 10.25 GHz. (b) One-photon Rabi oscillations of Psw 
with exponentially damped oscillation fits. Both the qubit 
Larmor frequency /qb and the microwave frequency /mwi 
are 10.25 GHz. The external flux is >I>qb/$o = 1.4944. (c) 
Four-photon Rabi oscillations when /qb = 41.0 GHz, /mwi ~ 
10.25 GHz, and $qb/$o = 1.4769. (d) The microwave am- 
plitude dependence of the Rabi frequencies nRabi/27r up to 
four-photon Rabi oscillations. The solid curves represent the- 
oretical fits. 



otherwise there should be no output voltage pulse. By 
repeating the measurement 8000 times, we obtain the 
SQUID switching probability Psw, which is directly re- 
lated to Pe{tp) for the dc readout pulse with a proper 
amplitude. For (f>qb > 1.5$07 -Psw is directly related to 
l-Po(ip). 

We performed a spectroscopy measurement of the 
qubit with long (50 ns) single-frequency microwave 
pulses. We observed multi-photon resonant peaks ($qb < 
1.5$o) and dips (<i>qb > l-5<i>o) in the dependence of Psw 
on /mwi at a fixed magnetic flux $qb. We obtained 
the qubit energy diagram by plotting their positions as 



a function of $qb/$o (Fig- Ela)). We took the data 
around the degeneracy point 'I'qb « l.S^o by applying 
an additional dc pulse to the microwave line to shift ^qb 
away from 1.5<I>o just before the readout, because the dc- 
SQUID could not distinguish the qubit states around the 
degeneracy point. The top solid curve in Fig.|2Ia) repre- 
sents a numerical fit to the resonant frequencies of one- 
photon absorption. From this fit, we obtain the qubit 
parameters Ej/h = 213 GHz, A/2tt = 1.73 GHz, and 
a — 0.8. The other curves in Fig. EJa) are drawn by 
using these parameters for ni = 2, 3, and 4. 

Next, we used short single-frequency microwave pulses 
with a frequency of 10.25 GHz to observe the coher- 
ent quantum dynamics of the qubit. Figures Efb) 
and (c) show one- and four-photon Rabi oscillations 
observed at the operating points indicated by arrows 
in Fig. Ufa) with various microwave amplitudes Vmwi- 
These data can be fitted by damped oscillations oc 
exp(— tp/Td) cos(5lRabiip), except for the upper two 
curves in Fig. |21[b). Here, ip and are the microwave 
pulse length and qubit decay time, respectively. To ob- 
tain Ofiabij we performed a fast Fourier transform (FFT) 
on the curves that we could not fit by damped oscilla- 
tions. Although we controlled the qubit environment, 
there were some unexpected resonators coupled to the 
qubit, which could be excited by the strong microwave 
driving or by the Rabi oscillations of the qubit. We con- 
sider that these resonators degraded the Rabi oscillations 
in the higher Vmwi range of Fig.EJb). FigureE^d) shows 
the Vmwi dependences of 17Rabi/27r up to four-photon 
Rabi oscillations, which are well reproduced by Eq. ||2Jl. 
Here, we used only one scaling parameter a(10.25 GHz) 
— 0.013 defined as a(/Mwi) = 4giQ!i/Li;MWi Vmwij be- 
cause it is hard to measure the real amplitude of the mi- 
crowave applied to the qubit at the sample position. The 
scaling parameter a (/mwi) reflects the way in which the 
applied microwave is attenuated during its transmission 
to the qubit and the efficiency of the coupling between 
the qubit and the on-chip microwave line. In this way, 
we can estimate the real microwave amplitude and the 
interaction energy between the qubit and the microwave 
2%iQ!i by fitting the dependence of r2Rabi/27r on Vmwi- 
These results show that we can reach a driving regime 
that is so strong that the interaction energy 2hgiai is 
larger than the qubit transition energy hiOq^. 

We have also performed experiments with two mi- 
crowave friequencies /mwi and /mw2- First, we car- 
ried out a spectroscopy measurement by using long (50 
ns) two-frequency microwave pulses. In addition to reso- 
nances caused by the multi-photon absorption processes 
at multiples of each microwave frequency (/qb = j^i/mwij 
?^2/mw2), we also clearly observed those due to paramet- 
ric processes (/qb = |/mwi ± /mw2|) (not shown). 

We next investigated the coherent oscillations of the 
qubit through the parametric processes by using short 
two- frequency microwave pulses. Figure |3Ia) [(b)] shows 
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FIG. 3: Experimental results with two-frequency microwave 
pulses, (a) [(b)] Two-photon Rabi oscillations due to a para- 
metric process when /qb = /mw2 + [— ] /mwi . The solid curves 
are fits by exponentially damped oscillations, (c) [(d)] Rabi 
frequencies as a function of Vmwi , which are obtained from 
the data in Fig. |3[a) [(b)]. The dots represent experimental 
data when Vmw2 = 16.9, 23.5, 33.0, and 52.0 [50.1, 62.9, 79.1, 
and 124.7] mV from the bottom set of dots to the top one. 
The solid curves represent Eq. I^. The inset is a schematic 
of the parametric process that causes two-photon Rabi oscil- 
lation when /qb = /mw2 + [-] /mwl 

the Rabi oscillations of Pgw when the qubit Larmor fre- 
quency /qb = 26.45 [7.4] GHz corresponds to the sum 
of the two microwave frequencies /mwi — 16.2 GHz, 
/mw2 10.25 GHz [the difference between /mwi = 
11.1 GHz and /mw2 = 18.5 GHz] and the microwave 
ampHtude of MW2 Vmw2 was fixed at 33.0 [50.1] mV. 
They are well fitted by exponentially damped oscilla- 
tions oc exp(— tp/Td) cos(f2Rabiip)- The Rabi frequen- 
cies obtained from the data in Fig. 3(a) [(b)] are well 
reproduced by Eq. (PJ without any fitting parameters 
(Fig. I^Jc) [(d)]). Here, we used A, which was ob- 
tained from the spectroscopy measurement (Fig. |2Ja)) 
and used a(10.25 GHz) = 0.013 and a(16.2 GHz) = 
0.0074 [a(ll.l GHz) = 0.013 and a(18.5 GHz) = 0.0082], 
which had been obtained from Rabi oscillations by using 
single-frequency microwave pulses with each frequency. 
Those results provide strong evidence that we can achieve 
parametric control of the qubit with two-frequency mi- 
crowave pulses. 

In summary, we investigated nonlinear responses in the 



superconducting flux qubit. First, we observed multi- 
photon Rabi oscillations caused by up to four-photon 
transitions by using single-frequency microwave pulses. 
The microwave amplitude dependences of the Rabi fre- 
quencies are well reproduced by Bessel functions derived 
from a semi-classical model. Furthermore, we success- 
fully demonstrated parametric control of the qubit by us- 
ing two-frequency microwave pulses. We observed Rabi 
oscillations of the qubit caused by parametric transi- 
tions when /qb = I /mwi i /mw2|- The Rabi frequencies 
as a function of the microwave amplitudes are well de- 
scribed by the product of two Bessel functions. These 
results indicate that the flux qubit offers a good testing 
ground for exploring quantum nonlinear phenomena in a 
macroscopic quantum object. Furthermore, these multi- 
photon processes observed in our experiment widen the 
frequency range of microwaves for controlling flux qubit. 

We thank T. Kutsuzawa, F. Deppe for useful discus- 
sions and for helping with the experimental setup. We 
also acknowledge useful discussions with J. Johansson, H. 
Nakano, M. Thorwart, Y. Yamamoto, J. E. Mooij, C. J. 
P. M. Harmans, and Y. Nakamura. This work has been 
supported by the CREST project of the Japan Science 
and Technology Agency (JST). 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computa- 
tion and Quantum information (Cambridge Univ. Press, 
Cambridge, 2000). 
[2] Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Nature 398, 

786 (1999). 
[3] D. Vion et al, Science 296, 886 (2002). 
[4] Y. Yu, S. Han, X. Chu, S.-I. Chu, and Z. Wang, Science 

296, 889 (2002). 
[5] J. M. Martinis, S. Nam, J. Aumentado, and C. Urbina, 

Phys. Rev. Lett. 89, 117901 (2002). 
[6] J. E. Mooij et al, Science 285, 1036 (1999). 
[7] I. Chiorescu, Y. Nakamura, C. J. P. M. Harmans, and J. 

E. Mooij, Science 299, 1869 (2003). 
[8] T. Kutsuzawa et al., Appl. Phys. Lett, (to be published); 

cond-mat/0501592 

[9] B. L. T. Ploude et a/., |cond-mat/0501679| 
[10] L. M. K. Vandersypen at., Nature 414, 883 (2001). 
[11] G. Yusa, K. Muraki, K. Hashimoto, and Y. Hirayama, 

Nature 434, 1001 (2005). 
[12] T. Yamamoto, Y. A. Pashkin, O. Astafiev, Y. Nakamura, 

and J. S. Tsai, Nature 425, 941 (2003). 
[13] R. McDermott et al, Science 307, 1299 (2005). 
[14] A. Wallraff et al, Nature 431, 162 (2004). 
[15] Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Phys. Rev. 

Lett. 87, 246601 (2001). 
[16] A. WaUraff, T. Duty, A. Lukashenko, and A. V. Ustinov, 

Phys. Rev. Lett. 90, 037003 (2003). 
[17] S. Saito et al, Phys. Rev. Lett. 93, 037001 (2004). 
[18] C. A. Sackett et al. Nature 404, 256 (2000). 
[19] C. Cohen- Tannoudji, J. Dupont-Roc, and G. Grynberg, 
Atom-Photon Interactions (Wiley, New York, 1992). 



